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Abstract 

We examine the LHC discovery potential for new (5 = 2/3 quark singlets T in 

the process gg,qq — > TT W~^bW~b, with one W boson decaying hadronically 

and the other one leptonically. A particle-level simulation of this signal and its 

I main backgrounds is performed, showing that heavy quarks with masses of 500 

CN . 

GeV or lighter can be discovered at the 5a level after a few months of running, 
^ , when an integrated luminosity of 3 fb is collected. With a luminosity of 100 

, fb~^, this process can signal the presence of heavy quarks with masses up to 

^ , approximately 1 TeV. Finally, we discuss the complementarity among TT, Tj 

production and indirect constraints from precise electroweak data in order to 
^ . discover a new quark or set bounds on its mass. 

i> : 

oo ■ 

^ '■ 1 Introduction 
vn ■ 

■ The Large Hadron Collider (LHC) will be a powerful tool to explore energies up to 

PhI the scale of a few TeV. It is expected to provide some striking evidence of new physics, 

Qh! for instance of a light Higgs boson, in its first months of operation [1,2]. Among 

^ , many promising possibilities for the discovery of new particles, LHC will offer an ideal 

> ■ environment for the production of heavy quarks. New quarks of either charge can 

^ ■ be copiously produced in pairs through QCD interactions, namely via gluon fusion 

' and quark-antiquark annihilation, if there is available phase space [3,4]. Up-type 

quarks T can also be produced in association with light jets, e.g. in the processes 
qb — » g'T, q'b qT (here and throughout this Letter q = u,c, q' = d,s), provided 
their mixing with the bottom quark is sizeable. New interactions may also bring about 
further production mechanisms. The prospects for heavy quark detection depend on 
the production processes (with their respective cross sections) as well as on the decay 
modes (and their relevant backgrounds), which are distinctive of the Standard Model 
(SM) extension considered. 
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The presence of a fourth sequential generation is disfavoured by naturalness argu- 
ments^ and precision electroweak data, which leave a small window for the new quark 
masses consistent with the experimental measurement of the S, T, U parameters [5]. 
On the other hand, heavy SU(2)l quark singlets with charges Q = 2/3 or Q = —1/3 
can exist with a moderate mixing of order 10^^ — 10^^ with the SM quarks. Here we 
are concerned with the first possibility. Models with large extra dimensions with for 
example tn in the bulk predict the existence of a tower of Q = 2/3 singlets T^"^. If 
there is multilocalisation the lightest one, T^li can have a mass of 300 GeV or larger, 
and a sizeable mixing with the top quark [6]. (The class of extra-dimensional models 
having a light Q = 2/3 singlet mixing with the top quark is enlarged when corrections 
localised on the branes to the kinetic terms of fermions and bosons are taken into 
account [7].) Little Higgs models [8] include in their additional spectrum an up-type 
singlet, which is expected to have a mass of 1 TeV or larger. Quark singlets also appear 
in some grand unified theories [4,9]. Their effects in low energy and top physics have 
already been studied [10]. In this Letter we address their direct observation at LHC 
through pair production gg, qq TT [11]. 

We note that for heavy quark masses > 800 GeV and a coupling to the bottom 
quark Vxh of the size suggested by the experimental measurement of the T parameter, 
single T production pp Tj has a larger cross section than pair production and can 
then explore larger mass scales [12,13]. Therefore, Tj production will eventually set 
more stringent limits (albeit dependent on Vrb) on heavy quark masses if a positive 
signal is not observed. However, two important points have to be remarked: («) the 
Tj cross section is proportional to |Vr6p, hence for small mixings this process becomes 
less relevant; (n) pair production has the best sensitivity to the presence of new quarks 
having masses of several hundreds of GeV. If new quarks exist in this mass range, TT 
production would allow to observe a signal in a rather short time. 

In the following we briefiy review the mixing of the new quark, its interactions and 
decay modes. After summarising the relevant aspects of the signal and background 
generation, we will present our results for quark masses of 500 GeV and 1 TeV. Finally, 
the relation between TT, Tj production and indirect constraints from the T parameter 
will be discussed. 

^For a fourth quark generation, anomaly cancellation requires the simultaneous presence of a lepton 
doublet. LEP measurement of the Z invisible width sets the number of light neutrino species to three, 
and additional neutrinos must be heavier than 45 GeV [5], in sharp contrast with the smallness of the 
hght neutrino masses rrii, < 1 eV. 
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2 SM extensions with Q = 2/^ singlets 



The addition of two 811(2)^ singlet fields to the quark spectrum modifies the 
weak and scalar interactions involving Q = 2/3 quarks. (We denote weak eigenstates 
with a zero superscript, to distinguish them from mass eigenstates which do not bear 
superscripts.) Using standard notation, these interactions read 

Cw = —^[u-i^VPLdW+ + dYV^PLuW-], 
v2 



XPr - 



4 
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Ch = -^u[M''XPl + XM''Pr]uH, (1) 

where u = {u,c,t,T), d = {d,s,b) and Pr^l = (1 ±75)/2. The extended Cabibbo- 
Kobayashi-Maskawa (CKM) matrix V is of dimension 4x3, X — VV^ is a non-diagonal 
4x4 matrix and Al" is the 4x4 diagonal up-type quark mass matrix. The new mass 
eigenstate T is expected to couple mostly with third generation quarks t, b, because 
Tl, preferrably mix with t^, t%, respectively, due to the large top quark mass. Vrb 
is mainly constrained by the contribution of the new quark to the T parameter [10,14], 

T- .^^' 2 {\VTb\'[0+{yT,yb)-e+{yt,y,)] - \XtT\'e+{yT,yt)} , (2) 

where iV^ = 3 is the number of colours, ?/j = {ml/Mz)'^, ml being the MS mass of the 
quark i at the scale Mz, l^tTp — 1^x6^(1 - \VTh\^) and [14] 

0+{yuy2) = yi + y2 log — . (3) 

yi - vi z/2 

The experimental measurement T = 0.12 ± 0.10 [15], obtained setting U = 0, implies 

T < 0.28 with a 95% confidence level (CL), and the corresponding limit \VTb\ < 

0.26 - 0.18 for mr = 500 - 1000 GeV (see also Ref. [16]).^ Mixing of with <, 

c°, especially with the latter, is very constrained by parity violation experiments and 

the measurement of Rc and at LEP, respectively [17, 18], implying small X^t, 

XcT- The charged current couplings with d, s must be small as well, \VTd\-i \ Vts\ ~ 0.05, 

because otherwise the new quark would give large loop contributions to kaon and B 

physics observables [10]. Therefore, |l^d|,|l^s| <S l^rftl and I-'^ctI <S l-'^tTl- 

^The new quark contribution to U is much smaller, of order 10~^ [10], thus it makes sense using 
this value for T. If we take T = — 0.17± 0.12 [5] the limits obtained are much stronger, T < 0.027 and 
thus \ VTb\ < 0.08 — 0.06 for rriT = 500 — 1000 GeV. We will consider both possibilities in our analysis. 
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In specific models there may be additional interactions, e.g. mediated by new gauge 
bosons, giving further contributions to experimental observables. These extra terms 
might (partially) cancel the ones from the new quark, loosening the constraints on 
its couplings. Although new interactions may modify the allowed range for Vrb, it is 
unlikely that they alter the above hierarchy. Then, the relevant decays of the new 
quark are T — > Zt, Ht, with partial widths 
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approximately equals mT/2 for mr ^ m^, M. For a heavy T and a light Higgs, we have 
Br(r ^ W+h) ~ 0.5, Br(r ^ Zt) ~ 0.25, Br(r ^ Ht) ~ 0.25. To our knowledge, 
there are not experimental searches for new (5 = 2/3 quarks giving lower bounds for 
their masses. However, for the expected production cross sections and decay branching 
ratios, it seems that quarks with masses around 200 GeV ought to be visible with 
present Tevatron Run II data [19]. 

The decays T — > — > t^l~W'^h, £ — e, /i give a cleaner final state than T — > W'^b, 
but with a branching ratio 30 times smaller. It has already been found that the channel 
T — > W'^b (T — > W~b), with W — > ii/, gives the best discovery potential in single T 
production [12]. In TT production we select the final states TT — > W'^bW~b, with one 
W boson decaying leptonically and the other one hadronically. The larger cross section 
in this decay mode allows to obtain a better statistical significance for the signal, while 
the backgrounds can be greatly reduced with kinematical cuts. 
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3 Signal and background simulation 



The main backgrounds for the TT signal 

gg,qq^Tf ^W+bW-b^e+iybqq'b , i = e, fM (6) 

are given by tt, Wbbjj, Zbbjj and tbj production, 

99: qq^ti^ W+b W'b tub qq'b , 
pp — > Wbbjj ii^bbjj , 
pp Zbbjj £^£~bbjj , 

pp tbj W+bbj e+ubbj . (7) 

The charge conjugate processes are understood to be summed in all cases. In the 
background evaluation we do not consider final states with r leptons (which can decay 
leptonically r — ez/T-z/g, t /uz/^z/^) because the electrons and muons produced in r 
decays are softer, and in our analysis we eventually require e, ^ with high transverse 
momenta. Wjjjj and Zjjjj production are reduced to negligible levels with the 
requirement of two b tags, which suppresses their contributions by a factor ~ 10~^. The 
signal and the ti, tbj backgrounds are evaluated with our own Monte Carlo generators, 
including all finite width and spin effects. We calculate the matrix elements using 
HELAS [20] with running coupling constants evaluated at the scale of the heavy quark 
T or t. Wbbjj and Zbbjj are calculated with ALPGEM [21]. The bottom quark mass 
mj, — 4.8 GeV is kept in all cases, and we take Mh — 115 GeV. We use structure 
functions CTEQ5L [22], with ^ s for TT, tt and tbj, and = M^,z + Pt^^^z for 
Wbbjj, Zbbjj, being y/§ the partonic centre of mass energy, and pxwz ^^e transverse 
momentum of the gauge boson. ^ 

The events are passed through PYTHIA 6 . 228 [24] as external processes to perform 
hadronisation and include initial and final state radiation (ISR, FSR) and multiple 
interactions. We use the standard PYTHIA settings except for b fragmentation, in 
which we use the Peterson parameterisation with ei, — —0.0035 [25]. A fast detec- 
tor simulation ATLFAST 2.60 [26], with standard settings, is used for the modelling 
of the ATLAS detector. We reconstruct jets using a standard cone algorithm with 

^We find that TT cross sections are 16 — 18% larger (for tut = 500 — 1000 GeV) when evaluated 
with MRST 2004 structure functions [23] and their corresponding value of as{Mz)- Assuming that 
background cross sections in the kinematical region of interest (large transverse momenta and invariant 
masses) scale by the same rate, this would amount to a 8 — 9% increase in the statistical significance. 
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AR = ^{At])^ + (A0)2 = 0.4, where 7] is the pseudorapidity and the azimuthal 
angle. We do not apply trigger inefficiencies and assume a perfect charged lepton 
identification. The package ATLFASTB is used to recalibrate jet energies and perform b 
tagging, for which we select efficiencies of 60%, 50% for the low and high luminosity 
LHC phases, respectively. 

4 Numerical results 

The hadronised events are required to fulfill these two criteria: (a) the presence of one 
(and only one) isolated charged lepton, which must have transverse momentum pt > 20 
GeV and \ri\ < 2.5; (b) at least four jets with pt > 20 GeV, \ri\ < 2.5, with exactly two 
b tags. The cross sections times efficiency of the five processes after these pre-selection 
cuts are collected in Table [U using a 60% b tagging rate. The events are produced 
without kinematical cuts at the generator level in the case of TT, ti, while for the other 
processes we set some loose cuts, less restrictive than the ones used after hadronisation, 
which do not bias the calculation. For tbj we only require pseudorapidities |?7| < 3 for 
6, j. For Wbbjj we set pt > 15 GeV and |?7| < 3 for the charged lepton, the b quarks 
and the jets, and lego-plot separations ARjj, ARbj, AR^b > 0.4, ARib, ARij > 0.2. For 
Zbbjj we require pt > 15 GeV, |?7| < 3 for 6 quarks and jets, \ri\ < 10 for the charged 
leptons, and ARjj, AR^j, AR^b > 0.4. 



Process 


a X eS 


Tf (500) 


44.9 fb 


Tf (1000) 


0.638 fb 


ti 


18.8 pb 


Wbbjj 


1.23 pb 


Zbbjj 


246 fb 


tbj 


710 fb 



Table 1: Cross sections of the TT signal (with my = 500, 1000 GeV) and its back- 
grounds after pre-selection cuts. 

The SM backgrounds are much larger than the signal, but they concentrate in 
the low transverse momenta region. To reduce them, it is useful to examine their 
dependence on the transverse momenta of the charged lepton p^^^, the fastest jet pj,'"^^^ 
and the fastest jet Pt'^^^, as well as the missing transverse momentum "^t and the 
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total transverse energy Ht = ^^Pt + ■ The kinematical distributions of these 
variables are shown in Fig. [H We display a weighted sum of the Wbbjj and Zbbjj 
processes so as to reduce the number of histograms, while the other backgrounds are 
shown separately. 

The TT signal can be discovered by the presence of peaks in the invariant mass 
distributions corresponding to the two decaying quarks. In order to reconstruct their 
momenta we first identify the two jets ji, j2 from the W decaying hadronically. The 
first one ji is chosen to be the highest pt non-6 jet, and the second one j2 as the non-6 
jet having with ji an invariant mass closest to M\y. The missing transverse momentum 
is assigned to the undetected neutrino, and its longitudinal momentum and energy are 
found requiring that the invariant mass of the charged lepton and the neutrino is the 
W mass, {pi + PuY = ^w- This equation yields two possible solutions. In addition, 
there are two different pairings of the two b jets to the W bosons decaying hadronically 
and leptonically, giving four possibilities for the reconstruction of the heavy quark 
momenta. We select the one giving closest invariant masses for the quarks 

decaying hadronically and semileptonically. Their kinematical distributions are shown 
in Fig. [21 In our calculations we have set Vt6 = 0.2, 0.1 for rriT = 500, 1000 GeV, 
respectively, yielding the total widths F^^ = 2.80, 6.16 GeV . The cross sections are 
independent of Vxh and for Ft of these sizes the broadness of the mass distributions 
too. 

We point out that in our signal calculation we have not included other T production 
processes giving the same experimental signature of one charged lepton, four jets (with 
two b tags) plus missing energy. Such processes do not constitute a background (they 
are absent in the SM) but instead increase the signal cross section. One example is 
TT production in the decay channel TT — * ZtW^b uuW^bW'b, with one W 
boson decaying hadronically and the other one leptonically. This process has a cross 
section 10 times smaller than the one in Eq. Other possible TT decay channels 
are TT ZtW-b bbWHW-b, TT HtW-b ^ bbW+bW-b (assuming a light 
Higgs boson) with two b quarks mistagged. Their contributions represent a ~ 7% and 
~ 40% increase, respectively, in the total cross section. Nevertheless, neither of the 
three processes mentioned yields peaks in the m^'^, n%p' invariant mass distributions, 
as reconstructed here for the TT W^bW~b signal, and their contributions are not 
likely to be detectable due to the uncertainty in the SM background normalisation. 
The same comments apply to Tbj production with radiation of an extra hard jet. 
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Figure 1: Transverse momentum of: (a) the fastest jet; (b) the fastest b jet; (c) the 
charged lepton. Missing transverse momentum (d); total transverse energy (e). The 
histograms are normalised to a total number of 2000 events. 
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Figure 2: Reconstructed masses of the heavy quarks decaying hadronically (a) and 
semileptonically (b). The histograms are normahsed to a total number of 2000 events. 



4.1 Results for uit = 500 GeV 

The large background cross sections make it convenient to introduce further kinematical 
cuts at the generator level to reduce the number of events processed with PYTHIA and 
ATLFAST. We require the presence of a charged lepton with pt > 30 GeV, a jet with 
Pt > 200 GeV and, for Wbbjj and Zbbjj, one b quark with pt > 100 GeV. This last 
cut does not bias the sample because in these two processes the two non-6 jets mostly 
originate from light quarks and gluons, for which the b mistag probability is very low. 
Thus, the 6-tagged jets correspond to the b quarks most of the time. The kinematical 
cuts used to reduce backgrounds are 

^,,max ^ GeV , p?'"^''" > 150 GeV , 

Pt"" > 50 GeV , 50 GeV < A< 600 GeV , 

Ht > 1000 GeV . (8) 

The cut 600 GeV is useful because ti production with large invariant masses is 
sometimes associated to very large ^t, in contrast with the signal. We also note that 
with these requirements the charged lepton and the hardest jet provide a trigger in 
the low luminosity LHC phase. The cross sections at the generator level are listed 
in the first column of Table [21 mainly for informative purposes. The second column 
represents the number of events A^o = KaC simulated, taking a luminosity of 10 fb~^ 
and including the rescaling factors K as explained in the appendix. The figures in these 
two columns corresponding to different processes should not be compared, since they 
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are obtained with different initial cuts in tfie event generation. Instead, tfie number 
of events iVcut surviving the selection criteria in Eq. ^ reflect the relative size of the 
processes after cuts. They are shown in the third column. (The size of the signals and 
backgrounds before cuts can be read from Table [H) 



Process 


a 




iVcut 


^peak 


Tf (500) 


204 fb 


2700 


272 


173 


tt 


5590 fb 


70000 


1609 


240 


Wbbjj 


928 fb 


16000 


287 


65 


Zbbjj 


364 fb 


7200 


39 


10 


fbj 


626 fb 


8300 


70 


11 



Table 2: For each process: cross sections a including cuts at the generator level; number 
of events simulated Nq; number of events N^ut passing the selection criteria in Eq. ((HI); 
number of events iVpeak passing the selection cuts which are in the peak regions. 

These kinematical cuts allow to detect the presence of the new quark in the invariant 
mass distributions m^'^, rrtp' as can be observed in Fig. El The number of events in 
the peak regions 

340 GeV < m^^^ < 660 GeV , 340 GeV < m!^P < 660 GeV (9) 

are displayed in the fourth column of Table El They give a statistical significance 
S/\/B = 9.6. A 5(T significance, needed to claim discovery, can be achieved with a 
luminosity C ~ 2.7 fb~^. These numbers only consider statistical uncertainties, assum- 
ing that the SM background can be normalised with the cross section measurements 
ouside the peak region in Eq. (jH)). Additionally, the trigger and charged lepton detec- 
tion efficiencies must be taken into account, what reduces the statistical significance 
by a factor ~ 0.95. 

4.2 Results for mr = 1 TeV 

We repeat the same analysis for a heavy quark with uit = 1 TeV, in this case choosing 
a b tagging rate of 50% at the high luminosity phase. The generator cuts are raised to 
Pt > 150 GeV for the charged lepton, pt > 250 GeV for the hardest jet and pt > 150 
GeV for the hardest b quark, the latter cut only for Wbbjj and Zbbjj production. The 
parton-level cross sections for the five processes are listed in Table El together with the 
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(b) 



Figure 3: Reconstructed masses of the heavy quarks decaying hadronically (a) and 
semileptonically (b), after the selection cuts in Eq. (jH]). The dashed hues correspond 
to the SM predictions, while the full lines represent the SM plus a new 500 GeV quark. 



number of simulated events No, corresponding to an integrated luminosity of 300 fb ^. 
The selection criteria used to reduce backgrounds are 

Pt" > 200 GeV , 50 GeV < ^t< 400 GeV , 

ift> 1800 GeV. (10) 

With these cuts the charged lepton and hardest jet provide a trigger for the event. The 
peak regions in this case are defined as 

800 GeV < 77^^^ < 1200 GeV , 800 GeV < m^^ < 1200 GeV . (11) 

The invariant mass distributions m^^'^, after the cuts in Eq. (fTIIll are shown in 
Fig. m The excess of events in the peak regions amounts to 9.4 standard deviations of 
the expected SM background. A 5(7 significance would be achieved with an integrated 
luminosity £ ^ 85 fb^^. With a simple rescaling it can be estimated that masses up 
to rriT = 1.1 TeV can be discovered with ba significance for C = 300 fb^^ and, if no 
signal is found, the 95% CL limit mx > 1-3 TeV can be set. 



5 Summary and discussion 

Up-type quark singlets with charge Q = 2/3 are predicted in some SM extensions, 
with masses ranging from few hundreds of GeV to several TeV. Their observation 
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Process 


a 




iVcut 




Tf (1000) 


2.89 fb 


1330 


70 


48 


tt 


778 fb 


294000 


208 


10 


Wbbjj 


66.8 fb 


34000 


132 


15 


Zbbjj 


48.0 fb 


28500 


19 


1 


tbj 


44.1 fb 


17500 


3 






Table 3: For each process: cross sections a including cuts at the generator level; number 
of events simulated Nq] number of events Nc^t passing the selection criteria in Eq. (fTTH) : 
number of events A^peak passing the selection cuts which are in the peak regions. 
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Figure 4: Reconstructed masses of the heavy quarks decaying hadronically (a) and 
semileptonically (b), after the selection cuts in Eq. (|TTHl . The dashed lines correspond 
to the SM predictions, while the full lines represent the SM plus a new 1 TeV quark. 



would represent not only a clear new physics signal but also an important confirmation 
for these models. We have shown that for rriT = 500 GeV a 5a statistical significance 
would be attained already with 3 fb~^ of integrated luminosity, which can be collected 
after few months of LHC operation in its first phase. For new quarks in this mass 
range, TT production provides the best signal of their presence, allowing a prompt 
discovery if they exist. With an integrated luminosity of 300 fb~\ TT production may 
discover a new quark with < 1.1 TeV, or set a 95% CL bound rriT > 1-3 TeV, 
independent of Vrb, if no signal is observed. 

We point out that if a fourth quark generation {T,B) exists with rriT < ms the 
dominant decay of the up-type quark is T ^ Wb, giving the same signal studied 
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here. (T — * Zt, Ht are forbiden at the tree level by the vanishing of the flavour- 
changing neutral coupling Xit-) The results obtained for an up-type singlet can then 
be straightforwardly applied to a fourth generation quark, multiplying the statistical 
significances in section |31 by a factor of four. If ttit < ttib, a fourth generation quark T 
with rriT < 1-3 TeV could be discovered with 300 fb~^ of integrated luminosity, and the 
95% CL bound rriT > 1-5 TeV could be set if they are not observed. If i? is lighter than 
T, the decay T — » W^B^ with W on its mass shell if my > rns + M^, is open. Hence, 
the branching ratio of T — > Wb depends on Vrb, and model-independent predictions 
cannot be made. We also remark that in case that a new quark T is discovered without 
a Q = — 1/3 partner, the experimental search for the decays T — > Zt, Ht can determine 
if T is a S'[J{2)l singlet or belongs to a doublet (in which case its partner B should be 
heavier and undetected). 

New (5 = 2/3 quark singlets can be produced in association with light jets as well, 
mainly in the processes ub dT, db — ^ uT. The cross sections for Tj production are 
quadratic in \Vtb\, but on the other hand they do not decrease with my as quickly as 
for TT. We plot in Fig. El (a) the cross sections for TT and Tj production for different 
heavy quark masses. For Tj we select a fixed coupling Vrt = 0.1 as well as an rriT- 
dependent coupling suggested by the experimental central value T = 0.12 (obtained 
for U = 0) and Eq. (|2|. The Vrb values derived from T = 0.12 are shown in Fig. [ni(b). 
Of course, in models beyond the SM there may be additional contributions to oblique 
corrections, thus we take the experimental measurement of T only as a hint on the size 
of Vrb- In Fig. El (b) we also plot the Little Higgs model relation Vrb = rnt/rriT and the 
two 95% CL upper limits on |Vt6| obtained from 

T < 0.28 (95% CL , U = 0) , 

T < 0.027 (95% CL , U arbitrary) , (12) 

respectively. We note that that the assumption Vrb = mt/rriT potentially conflicts with 
the T parameter measurement for rriT < 900 GeV, even using the less restrictive bound 
T < 0.28. 

The discovery potential of Tj production can be estimated from existing analyses. 
This process, with T Wb iub, gives a 21.5cr signiflcance for my = 1 TeV and 
Vrb = rrit/mT — 0.175 [13]. We make the reasonable assumption that for different 
T masses the signal to background ratio S/B in the kinematical region of interest 
(with high pt and invariant masses ~ m-j-) remains approximately constant (obviously 
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Figure 5: (a) Cross sections for TT production (full line) and Tj production, in the 
latter case for Vrt = 0.1 (dotted line), and for Vrb derived from the T parameter 
(dashed line), (b) Upper bounds on |VTf,| and values suggested by the T parameter 
(black line) and the Little Higgs relation Vrt = mt/rriT (grey line). 



keeping equal CKM factors iV^feP for the signal).'' Requiring a statistical significance 
S/Vb = 5 sets a lower limit on the coupling Vrb for each value. These limits 
are plotted in Fig. IHl (a), together with the discovery limit < 1.1 TeV for T pair 
production. We also include the 95% CL bounds from the T parameter in Eq. (fT2|l . 
We point out that if T < 0.027 is enforced the discovery reach of TT production is 
higher than for Tj, since for > 700 GeV the Vrt values needed for 5cr discovery 
in Tj production are not allowed. Assuming the less restrictive limit T < 0.28, there 
is a region (light shaded area in the figure) where the new quark can be discovered in 
single T but not in T pair production. Conversely, in the dark shaded area the new 
quark can be discovered in TT production but not in single T processes. 

The discovery of a new Q = 2/3 quark singlet would certainly be a rather important 
achievement towards the understanding of the fiavour structure of the SM, and might 
help explain the largeness of the top quark mass [8,27]. On the other hand, the non- 
observation of new quarks at LHC would also be interesting on its own. In such case, 
the combined bounds obtained from single and T pair production and the T parameter 
would restrict m^, |Vrb| to lie in the shaded area in Fig. IHl (b). (If we use the more 

^This assumption is justified by the decrease of the tails in the transverse momenta and invariant 
mass distributions of the SM background. More optimistic extrapolations of the SM background cross 
section, e.g. assuming that it decreases faster than the Tj signal, would lead to higher discovery limits 
on T masses, as the ones given in Ref. [12]. 
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Figure 6: (a) ha discovery limits for the new quark (full lines), and indirect bounds 
from the T parameter, explained in the text, (b) Combined 95% CL bounds on mr, 
|Vt6| (shaded area) from TT, Tj production and the T parameter, if Q = 2/3 quark 
singlets are not observed at LHC. 

restrictive bound in Eq. (fT^ the allowed region is somewhat smaller, as can be seen in 
the figure.) In this area we have rriT > 1.3 TeV and \VTb\ < 0.13, the latter implying 
\Vtb\ > 0.991. For > 600 GeV, the couplings Vj-d, Vts are already very constrained 
by kaon and B physics measurements [10]. Therefore, the non-observation of a new 
quark would significantly improve the indirect limits on CKM matrix elements Vtd, Vts, 
Vtb within this class of models. 
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A Appendix: signal and background normalisation 

In addition to the processes listed in Eqs. ©,(111) there are other higher order processes 
contributing to the signal and backgrounds, namely final states including extra jets 
from QCD radiation. Indeed, tbj production only yields three jets at the partonic 
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level, and the fourth one required to pass our pre-selection criteria must be originated 
by radiation. These higher-order processes are approximately accounted for by PYTHIA 
showering, which generates hard extra jets by FSR. For example, in Wbbjj production 
some fraction of the events, when passed through PYTHIA showering, are converted 
into Wbbjj j or Wbbjjjj events, the additional jets with a high (> 50 GeV) transverse 
momentum. 

One possible method to take higher order processes into account is to generate 
them at the parton level, forbidding PYTHIA to radiate hard extra jets to avoid double 
counting [28] but allowing the soft and collinear ones. Instead, our approach is to 
normalise the numbers of events simulated A^'o by approximate correction factors K so 
that these figures correspond to the processes in Eqs. ((El),© plus higher order ones. 
The number of events simulated for a given process is then A^o = KaC, being a the 
cross section of the process and C the luminosity. The correction factor K is calculated 
in the Wbbjj example as follows: 

1. We generate Wbbjj events requiring high transverse momenta > 100 GeV for 
jets at the generator level, and a large separation Ai? > 0.6 among all partons. 

2. These events are passed twice through PYTHIA and ATLFAST, including FSR and 
without including it. In both cases ISR, multiple interactions and energy smear- 
ing are turned off, because we want to isolate the effect of FSR. 

3. We examine the number of events with four-jets (corresponding to the four initial 
partons) with pt > 80 GeV in both samples. Let us call these numbers n^, nf^^, 
respectively. K is then defined as the ratio between them, K = n1/nl^^. 

With this rescaling factor definition, the number of four-jet events after FSR corre- 
sponds to the one present at the parton level. The value of K depends on the cut 
values for the event generation (100 GeV) and the jet counting (80 GeV), and there- 
fore this procedure is approximate. However, for our purposes this simple and fast 
X-factor prescription seems sufficient. For the TT signal we obtain K{TT) ~ 1.3, 
K{TT) ^1.5 for heavy masses of 500 GeV and 1 TeV, respectively, and for the back- 
grounds K{tt) ~ 1.3, K{tbj) ~ 1.3, K{Wbbjj) ~ 1.7, K{Zbbjj) ~ 2.0. 
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